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We calculate exactly matrix elements between states that are not eigenstates of the quantum 
XY model for general anisotropy. Such quantities therefore describe non equilibrium properties of 
the system; the Hamiltonian does not contain any time dependence. These matrix elements are 
expressed as a sum of Pfaffians. For single particle excitations on the ground state the Pfaffians in 
the sum simplify to determinants. 



Spin systems are paradigmatic models for describing many phenomena in contemporary physics [1]. Often their 
main properties can be captured qualitatively resorting on approximated techniques. However, in many cases more 
refined approaches are required to obtain reliable results. A prime example are systems near or at a phase transition, 
where quantum fluctuations inhibit many standard routes from working (as e.g mean-field theory but also perturbation 
theory). At the critical point, in change, the system can remarkably be simplified by a then present large class of 
symmetries. Conformal field theory employs systematically this important property of the system at criticality and 
the corresponding dynamics can be integrated exactly in 1+1 dimensions. Much more difficulties arise when the 
system is far from either, criticality and accessibility to mean-field or perturbation theory. Fortunately there are 
many non-trivial systems for which the symmetry is large enough to allow the dynamics being integrated exactly 
even for generic values of the relevant couplings. Then, also the physics of the cross over from non-critical to critical 
regimes is accessible. Complete integrability constitutes the crucial property that even exact correlation functions are 
available. Important steps forward to this goal have become possible by the Quantum Inverse Scattering approach [2], 
more recently refined by Kitanine, Maillet, Slavnov and Terras [3] and Korepin and Gohman [4]. 

The quantum anisotropic XY chains are a relevant example of completely integrable models. The model was solved 
exactly by Lieb, Schultz and Mattis [5], Pfeuty [6] for isotropic cases and by Barouch, McCoy, and Dresden [7] for 
generic anisotropy. Also the correlation functions were intensively studied and analytic expressions for their asymp- 
totics (in time and space variables) were obtained [7,8]. The correlation functions were calculated at equilibrium 
and for time-dependent magnetic field. We perform an exact calculation of correlations between states that are not 
eigenstates of the model, and that therefore describe non- equilibrium properties of the model; we remark that the 
Hamiltonian instead does not contain explicit time-dependence. Our motivations come from condensed matter where 
quantum XY chains are particularly studied, even more intensively since recent interest in the phenomenon of deco- 
herence in suitably designed physical systems [9]; this latter kind of analysis is due, in turn, to the burst of interest 
in quantum information theory [10]. Such cross-over of interests originated a line of research investigating the inter- 
connection between condensed matter and quantum information. In particular it is intriguing to investigate whether 
it is possible to better characterize condensed matter states by looking at e.g. quantum correlations entanglement 
properties of their wavefunction. Already a number of interesting results in this direction have been obtained [11]- 



The present paper is laid out as follows. In the next section we present the models we discuss and review the exact 
solution from Refs. [5,7], already preparing relevant building blocks for computing off-equilibrium correlations. In 
section III we present known results connecting vacuum expectation values in fermionic theories with a generalized 
determinant structure, called the Pfaffian and their application to equilibrium correlation functions presented in [6,7]. 
Section IV contains the main result for non-equilibrium correlations and matrix elements of the presented models. 
After all we draw our conclusions. 



The system under consideration is a spin-1/2 ferromagnetic chain with an exchange coupling A in a transverse 
magnetic field of strength h. The Hamiltonian is H — hH s with the dimensionless Hamilton operator H s being 
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where S a are the spin-1/2 matrices (a = x, y, z) and N is the number of sites. We assume periodic boundary conditions. 
The anisotropy parameter 7 connects the quantum Ising model for 7=1 with the isotropic XY model for 7 = 0. In 
the interval < 7 < 1 the model belongs to the Ising universality class and for N = 00 it undergoes a quantum phase 
transition at the critical coupling A c = 1. The order parameter is the magnetization in x-direction, {S x ), which is 
different from zero for A > 1 and vanishes at and below the transition. On the contrary the magnetization along the 
z-direction, (S z ), is different from zero for any value of A. 

This class of models was diagonalized by means of the Jordan- Wigner transformation [5-7] that maps spins to one 
dimensional spin- less fermions with creation and annihilation operators c\ and c v It proved convenient to use the 
operators Ai = cj + c ( , Bi = c\ — c ; , which fulfill the anti-commutation rules 

{A u A m } = -{Bi,B m } = 2S lm , 

{A u B m } = 0. (2) 

In terms of these operators the Jordan- Wigner transformation reads 

1-1 

S? = -A x [] A S B S 
z s =i 

S y l ^-\B l l ^A 8 B s 

s=l 

Sf = - 1 -A l B l . (3) 

The Hamiltonian defined in Eq.(l) is bilinear in the fcrmionic degrees of freedom and therefore can be diagonalized 
by means of the transformation 

^ = 77^ 12 elkl ( akc i + z ^ c ( 4 ) 



with coefficients 

oik 



A/t — (1 + Acosfc) 
v /2[A^-(l + AcosA:)A fe ] 



o = 7Asinfc 

^^-(l + Acos^] ■ 

The Hamiltonian thereafter assumes the form 

H = i2^vhk-lJ2 Ak ( 6 ) 
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and the associated energy spectrum is 



A fe = + A cos kf + A 2 7 2 sin 2 k . 



Now, in order to calculate correlations out of equilibrium, we need to know the time dependence of the relevant 
operators. From Eq.(4) we obtain the spin-less fermion creation and annihilation operators in the Heisenbcrg picture. 
We have ^(t) = exp (— i Akt)r)l(0) and hence, using Eq.(4) and its inverse 

c j (t) = 12[^-j(t)c l -b l - j (t)cj} 



where the new coefficients are 



a x (t) = -U^cosfca: (e iAfct - 2^ sin A fe t) (7) 

* k 

2i x - 

b x (t) — —= 2 ^ sin kxakflk sinA fc t . (8) 
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In the limit 7 = the previous expressions simplify considerably. In this case the magnetization, i.e. the z- 
component of the total spin S z = 5^ • S?, is a conserved quantity. In terms of fermions this corresponds to the 

conservation of the total number of particles, N — J2j n j = J2j c ] c j- For 7 — ► and |A| < 1 we find that a.k — ► 
and (3k — > signfc. The energy spectrum is At = |1 + Acosfc| and the eigenstates are plane waves (the Hamiltonian 
corresponds to a tight binding model) 

c i(*) = ^ EE cosfc(/- jK^'q (9) 



(10) 



In this work we discuss vacuum expectation values and correlations in excitations of them. It is worthwhile noticing 
that different strategies are applied, depending on whether the vacuum is the ground state or the state with no 
particles (which we call the c- vacuum). Since it is cumbersome to calculate the time dependence of the vacuum itself, 
it is convenient to write the operators Ai and B\ in the Heisenberg picture. For the ground state instead the time 
dependence is trivial and the operators are taken in the Schrodinger picture. For both approaches we express the 
operators Ai and Bi in terms of the operators i]k and ry k 



^ = 7^E[^ + ^ 

v q 



-i ql 



-i ql 
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These are sufficient for the correlations in the ground state. For the calculation for the c-vacuum it proves to be 
convenient defining the following (redundant) Fourier transforms containing a*,, (3k, an d their combination Zk := 

a k + i/3fe 



a, 

25* 
23* (*) 



q 

l^d^V** 
q 



In these quantities we have 



where 



M*) :=a5*(t)-i»3xW • 

^(*) = E(c5/*j-i(*) + c^_,(t)) 

j 

^(*) = E( c ]w-i(i)-c j ft%(t)) • 

j 

/U x (i) = y ^2 d<je~ lA<!t (a^ cosgx + a q (3 q sinqx) . 
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III. CORRELATION FUNCTIONS AS PFAFFIANS 



It is known since 1952 that "vacuum" expectation values of a product of 2R fermionic fields 

(0| *i ■ ■ ■ * 2 fl |0> (22) 

can be written as a Pfaffian [20] . The entries of the Pfaffian structure are the contractions of two field operators 

(01*^10)=^, (23) 

The field operators \1/ are linear functional of fermionic creation and annihilation operators, where the "vacuum" |0) 
is that state annihilated by the annihilation operators. The Pfaffian is a type of generalized determinant form [20]. 
It is written in a triangular structure as 

(^) 7r -P7r(l) : 7r(2)-f7r(3),7r(4) • • • ^ > ir(2n-l),7r(2n) = 



I -Pi, 2 Pi, 3 ■ ■ ■ P\,R Pl,R+l Pl.R+2 
P2,3 ■ ■ ■ P2,R P2,R+1 P2.R+2 

Pr.r+i Pr,r+2 

Pr+1,R+2 

P2R-1,2R 

where S^n denotes all elements ir of the symmetric group £2™ which give ordered pairs; i.e. ir(2l — 1) < 7r(2Z)) and 
tt(21 — 1) < 7r(2m — 1) for I < m. We particularly make use of the known property that a Pfaffian can be expanded 
along "rows" or "columns", where the r-th row or column corresponds to all Pij with i = r or j = r. In analogy 
to matrix minors, we will call the minor Pfaffian Pi j = Pjj the Pfaffian of the above structure (24) when having 
canceled the i-th and j-th row. In terms of these minors the expansion reads 

2fl-l 

V2R= EH i+r+lp r,%< (25) 

i=l 

where i, r means that the indices are to be written in increasing order. It is worth noting that the r-th part of this 
expansion reflects all possible contractions with the field operator \l/ r performed in Eq.(22). 

There are two cases which we will study in this work: |0) being (i) the ground state, denoted by \g) and (ii) the 
c- vacuum, denoted by |JJ.). 



P2.2R 



Pr,2R 
Pr+1,2R 



(24) 



A. Ground state 



At equilibrium [5-7], a crucial simplification is that (AiA m ) g = —(BiB m ) g = Si m . This reduces the Pfaffian to a 
Toplitz determinant 



(S?S? +R ) g = s(a,a)\0 ... 


^1,1 


"1,2 







U R,1 


"iJ-1,2 ■ 

^R.2 




■■ ^R-1,R 
■■ ^R.R 





= (_)*(*-l)/2 s(a>a ) 


"1,1 ■ 
^R,! ■ 


■ ■ ^l.R 
• ■ ^R,R 
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with [21] 



- {A l+ ^B l+v ) g 



(27) 
(28) 



(A[B m ) g = 3m-!- The correlation functions (SfSf +R ) g and (SfSf +R ) g identically vanish, since a complete row and 
column in the corresponding matrices vanishes, respectively. It is worth noting that due to the translational invariance 
of the state the determinant above is of Toplitz type. As a consequence the asymptotics of the correlation functions 
can be extracted explicitly [7] applying the Szego theorem. 

IV. CORRELATION FUNCTIONS OUT OF EQUILIBRIUM 

As already mentioned, time-dependent correlation functions were derived in Ref. [7]. There, the time-dependence 
was explicitly induced into the Hamiltonian (time dependent external magnetic field). In contrast, we compute 
matrix elements of operators at non-equilibrium, meaning that the initial and final state are not eigenstates of the 
Hamiltonian; the resulting quantities are then time dependent although the Hamiltonian is not. First, we consider 
matrix elements in the c- vacuum |-||) (for generic 7 this is not an cigcnstate of the Hamiltonian; only for 7 = it 
coincides with the ground state), and in excitations on it and on the ground state. 

A. Correlations in the c- Vacuum 

Using Eqs. (19,20), leads to the following contractions as building blocks for the Pfaffians 

(Ai(t)B m (t))^ =^2fj,*_ ll j im _ j = 



Sim - 4y ( 2a lPq COS «( m ~ + a <A I 1 ~ 2 @q) sin ?( m ~ 0) sin2 K l 



(29) 



(Ai(t)A m (t)}^ = y^M.*-;Mj-m = Sim - 2i ?jr a q f3 q sin q(m - I) sin2A 9 t 



(30) 



(Bi(t)B m (t)) 4 = ^2fj,*_ l ni- j = -Si m - 2iy ^a 9 /?,sinqr( m -Z)sin2A g i 



(31) 



We are now ready to write down the two-point spin correlation functions, applying the results from the previous 
section 



(S?St> +R h = s(a,f3y 



r af3 
1,2 



jap 
2 1,R-1 



T af3 



T< 



L R-2,R-1 



a/3 



F- 



ja/3 
J R-2 

ra/3 
J R-1 



a/3 



R-2 



a/3 
± R-l 

E al3 



s-ia.fi 
"1,2 



sia/3 
"_R-2,2 
fiaf) 
"-R-l, 2 



K. 



a/3 



H. 



a/3 
2,3 



f^a/3 



s-ta/3 
U R-2,R 
f~<a(3 
U R-1.R 

D< R 
Kf 



H 



a/3 

R-1,R 



(32) 
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where s(x,x) = s(y,y) = l/A(-) R ^ R+1 ^ 2 , 

1%, = (A l+ „(t)A l+ „(t))^ 



JXX JXX 



H^ v = (Bi +fl -i(t)Bi +v _i(t))j). 

TSXX TTXX 

1V u — n l,v 

= (A,+„(t)B, +v _i(t))^ 

jpxx r^xx 

JjJXX _ QXX 



T}XX flXX 



jyy — jyy 



Hvv v = (B l+li (t)B l+v (t)U 



K yy = H yy 



G™„ = (A l+lt - 1 (t)B, +v {t)h 
rryy = rpv 

A 1 M.l 



E yy = QVV^ 

D yy = 



and s{x, y) = s{y, x) = -i /4(_)«(«-i)/2 ; 



jxy 


= {A l+fl {t)A l+v {t))^ 


nxy 


= (A l+ ^t)B l+v (t)) 4 


j*y 




j?xy 

M 




ffxy 


= (B l+fl (t)B l+v (t))^ 


fixy 


rrxy 
n 0,l 


r>xy 

V 


^0,1/ 


K x v 


n \,v 



jyx 


= 


-i(t)Ai +v 


nyx 




-i(t)Bi +v 


jyx 


= 


jyx 


pyx 
n 


= 


jyx 

1 n,R+i 


Ejyx 


= 


1 R,R+l 


Dl x 




f~iyx 
u R,u 


Kv x 




f-iyx 
^R+l,v 






-i(t)Bi +v 



(36) 



We note that a Pfaffian P can be written (up to a sign) as a determinant of the corresponding antisymmetric matrix 
A of dimension 2R x 2R [22] by [pfP] 2 = det A. Since the c- vacuum is transactional invariant this determinant is 
again of Toplitz type. Therefore, also here the asymptotics of the correlation functions could be extracted explicitly 
along the lines depicted in [7] . 



B. Matrix elements for excitations of the vacuum 



We now want to concentrate on expectation values for states which are not the vacuum. So let and be linear 
functionals in the creation and annihilation operators and let us calculate 

(C\ *! ■ ■ ■ V 2 r \C) := (0| C*i ■ ■ ■ ^>2rC A |0) (37) 

Performing all possible contractions in (37), we obtain (0| CC'' f |0) (0| \&i • • • ^>2R |0) + all possible contractions where 
C and C'^ are contracted with a pair of field operators, say *j and \fr,-. Thus we have to calculate 

P itj := (C\ |0) (0| \C) - (C\ |0) (0| * 4 \C) , (38) 

where we take i < j in order to avoid double counting of contractions. The sign coming from transporting the operators 
C and C'^ to the left of \&j and the right of tyj respectively is (— ) t+ - 7 ' +1 . In the remaining vacuum expectation value 
the field operators "J/j and ^ are missing, which corresponds to canceling the rows i and j in the original Pfaffian 
(24). Consequently, this expectation value is the minor Pfaffian P^lj, and we obtain 

2R-1 2R. 

(C\ *! ■ ■ ■ * 2 R \C) := 7> 2fl + Yl E H <+j+1 4;% = 

1=1 j=»+l 



2ii-l / i-1 27? \ 

E EH <+j+1 °-%+ E 

We note that this expression is the sum over Pfaffian expansions (25). Indeed, each element of the first sum is the 
expansion of a Pfaffian along the i-th row, in which Pjj = and Pjj — ► Pij, hence 

2i?-l 

(C|*i---* 2 fl|C'):= E^L (39) 

i=0 

where we defined 

Vfl:=V2R 
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and 



Pl,2 ■■ 


. 


P\,i+1 ■ 


P\.2R, 









■ Pi-\aR 






Pi,i+1 


Pi,2R 








P2R~1,2R 



(40) 



Actually, we found the correlations (C\S"S^\C'} expressed as a sum of Pfaffians. The generalization to operators 
Ct that are multi-linear in the annihilation and creation operators can be related to a sum of multi-row expanded 
Pfaffians [23]. 

In what follows, we will come back to the cases of the ground state and the c- vacuum, discussed in the previous 
section. As mentioned, in order to be able to explicitly extract the asymptotics of the correlations, the initial and 
final state have to be transactional invariant. In the following, the translational invariance is explicitly broken. 



1. Single hole excitations on the ground state 



We choose the final and initial state to be (C\ = (g\ cj/\/?8o =: (j\ and \C) = c k \g) / V^o =: \k), which is normalized 
due to (g\ cjc fe \g) — ^k-j ■ Then we have to calculate 



(AiB m ){ k = (j\ A l \g) (g\ B m \k) - (j\ B m \g) (g\ A x \k) 



In this case wc find 



(AiB m ) ik 
(AiA rn ) g 



Q33 m _,-933fc- 
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jk 



«0 



(BiB rn ) 



jk 



So 



m-j 



The contraction of C with C is 



\CC\0) 



So 



33 



o 



(41) 
(42) 
(43) 

(44) 



We now discuss the correlation functions (S?S? +R ). The only non-zero contributions come from contractions of C 
and C with one operator of type A and one of type B (since only vacuum expectations of an equal number of A's 
and -B's are non-zero as discussed before). That means that here the sum of Pfaffians simplifies to the following sum 
of determinants 



<SfSf>=5><{ 



with 



2>W = (-1) J 



1,1 

G~ act 
U 



R,l 



with 



1,R 



i,R 



R,R 



jk 



G™ v — (Ai +tl Bi +u -i) g 
G™ v = {A l+ ^B l+v ) 



jk 
9 



(45) 



(46) 
(47) 
(48) 
(49) 
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Analogously, for the correlation functions (SfSf +R ) and (SfSf +R ) the only non-zero contributions come from con- 
tractions of C and C with two operators of type B and two of type A, respectively. This again simplifies the sum of 
Pfaffians to a sum of determinants 



2R ' 



where 



vf R = | 



G* 



xy 

1,1 • ' 


. 


s-ixy 


r ,xy 
■ ■ ■ ^1,R+1 


xy 

R-2,1 ■ ■ 

xy 

R-1,1 ■ ■ 
. 


. 
. 


f-ixy 

^R-2,i+l ■ 
f-ixy 

^iJ-l.i+l ' 



r-ixy 
■ ■ ■ ^R-1,R+1 
... 









frxy 


... 

frxy 



(50) 



and this simplifies to 



with 



(*) 

R 



i=l 



with 



V 



(0 



{-I) 1 



sixy 
U l,l 


• ■ ^l.i-l 


fixy 


r xy 
■ ^1,-R+l 





• • ^R-i,i-i 
.. 


r xy 

^R-1,%+1 ■ 
frxy 


n xy 
■ ^K-l.ii+l 
fjxy 


f-ixy 




a 




frxy __ 









(51) 



(52) 
(53) 



In an analogous way we find 



i=i 



with defined as in (51), but 



~~ (Ai+ v -iBi + ^) g 



(54) 
(55) 



2. Single particle excitations on the ground state 



Alternatively we consider (C\ = (g\ c^\f^r] =: (j\ and \C) = c\ \g) /\/2lo =: \k) as final and initial state, which arc 
again normalized according to (g\ c^c k \g) = %lk-j ■ In this case the possible contractions with the operators C and C 



are 
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(^B m ) J g = ^ 3 — (56) 



The contraction of C with C is here 



Oo 



(B/B ro ) s = — , (58) 



(0|CC'|0> = ^ (59) 



3. Single particle excitations on the c-vacuum 

We take the final and initial state to be (C\ = {0\cj —: {j\ and \C) = c£ =: \k). In this case the Pfaffian V\ R is 
given by Eqs. (32)-(36), where (following Eq. 40) in the i-th row 



are replaced with 



{AiB m )f = (j\ A l m (4| B m \k) - (j\ B m m (4| A l \k) 



and in the same manner (4| A t A m |4) — » (AiA m )f, (4| B ; B m |4> — » (BiB m )f. We find 

(Ai(t)B m (t))f = - (fj, m -jfj,%_i + Mj-;Mm-fe) (60) 

(Ai(t)A m (t))f = Mj-JMfe-m - Mj-m/4-J (61) 

(Bi(t)B m (t))f = Hm-jnt_ k - fii-jfi* m _ k . (62) 

With these results, all spin-correlation functions can be calculated as long as (C\ Sf(t = 0) \C) = (C\ Sf(t = 0) \C) = 0. 

In this case it will remain zero during the evolution. This is satisfied if the parity symmetry of the Hamiltonian is 
not broken by neither the initial nor the final state. 



V. CONCLUSIONS 



We calculated exactly spin-spin correlations out of equilibrium. For excitations on the "vacuum" , they can be written 
as a sum of Pfaffians. For excitations on the ground state these Pfaffians reduce to determinants (see Eqs. (39), (40)). 
The result for particle and hole excitations on the ground state and the c-vacuum were based on different approaches, 
writing the fermionic field operators in the Schrodinger and Heisenberg picture, respectively. Comparing with the 
known eigenstate correlations, we remark that here (S"S" +R ) cannot be reduced to R x R Toplitz determinants. For 
the vacuum-correlation functions the Pfaffians instead, they can be related to 2R x 2R Toplitz determinants. For 
correlation functions in excited states (of the vacuum) the translational invariance of the system is explicitly broken 
and then the determinants are not anymore of Toplitz type. This last issue constitutes a further difficulty of the 
problem to find the asymptotics of the correlations since the Szego theorem cannot be applied. 

We have used these results explicitly for studying the dynamics of correlations and quantum information theoretic 
quantities like the entanglement in specific states [19] but is also a key ingredient for the study of transport properties 
of the system. One possible application of the exact results we found here is to study the quantum phase transitions 
(characteristic of this class of models) out of equilibrium. 
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